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A thermodynamic approach to rapid phase transformations within a diffuse interface in a binary 
system is developed. Assuming an extended set of independent thermodynamic variables formed 
by the union of the classic set of slow variables and the space of fast variables, we introduce 
finiteness of the heat and solute diffusive propagation at the finite speed of the interface advancing. 
To describe the transformation within the diffuse interface, we use the phase-field model which 
allows us to follow the steep but smooth change of phases within the width of diffuse interface. 
The governing equations of the phase-field model are derived for the hyperbolic model, model with 
memory, and for a model of nonlinear evolution of transformation within the diffuse-interface. The 
consistency of the model is proved by the condition of positive entropy production and by the 
outcomes of the fluctuation-dissipation theorem. A comparison with the existing sharp-interface 
and diffuse- interface versions of the model is given. 
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I. INTRODUCTION 

Arising from phase transformations, a classic free- 
boundary problem introduces a model of phase interface 
with zero thickness. Within this problem, a sharp dis- 
continuity in properties or a jump of fluxes and ther- 
modynamic functions occurs across the interface. The 
sharp-interface model has been a successful description of 
many physical phenomena in various systems 0. How- 
ever, the sharp-interface model has difficulties in describ- 
ing situations when interfacial thickness becomes com- 
parable with the characteristic length of the considered 
phenomenon, and when a topology of the interface be- 
comes complicated or multiply connected. In order to 
avoid these difficulties in the sharp-interface model, an 
alternative model with a finite interfacial thickness was 
suggested for explaining phase transformations Q . 

Historically, the first formulation of basic principles of 
diffuse interfaces was given by Poisson, Maxwell, and 
Gibbs (3 who suggested to consider interface as a re- 
gion with finite thickness in which a steep but smooth 
transition of physical properties of phases occurs. Lord 
Rayleigh, van der Waals, and Korteweg |4| applied ther- 
modynamical principles to develop gradient theories for 
the interfaces with non-zero thickness. Through the past 
century, ideas of diffuse interface given by these authors 
0, 0] were refined and applied in many physical phenom- 
ena (see, e.g., overviews in Ref. y5|)- 

The diffuse-interface formalism has been widely ap- 
plied to phase transformations in condensed media. Bor- 
rowing the formalism of the Landau theory of phase 
transitions [||, the first introduction of the diffuse inter- 
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face into the theory of phase transformations was made 
by Landau and Khalatnikov who labelled the differ- 
ent phases by an additional order parameter to describe 
anomalous sound absorption of liquid helium. In its 
well-known form, the formal variational approach was es- 
tablished by Ginzburg and Landau for phase transitions 
from the normal to the superconducting phase On 
the basis of this approach, diffuse-interface models with 
order parameters have been developed by Halperin, Ho- 
henberg and Ma for the theory of critical phenomena [j| , 
and by Allen and Cahn for antiphase domain coarsening 

The diffuse- interface model has been also developed for 
the description of phase transformations of the first or- 
der, especially, for the solidification phenomenon. The 
diffuse-interface model of solidification incorporates an 
order parameter in the form of a phase-field variable . 
The phase-field <E> has a constant value in homogeneous 
phases, e.g. $ = — 1 for unstable liquid phase which is 
transforming into the solid phase with $ = +1. Between 
these phases in the interfacial region, the phase field, <E>, 
changes steeply but smoothly from —1 to +1. In numer- 
ical solutions, it allows one to avoid explicit tracking of 
the interface and locate the interface at $ = As 
a particular case, the phase-field model is reduced to the 
sharp interface limits and adopts the major models 
of sharp interface (such as Hele-Shaw type models, classi- 
cal or modified Stefan problem, etc.). The phase-field <& 
is considered as an order parameter which is introduced 
to describe the moving interfacial boundary between the 
initially unstable phase and the final phase. 

Several thermodynamically consistent phase-field 
models have been proposed [3 0, 0, 0] . These in- 
clude models for transformation in a pure system [l4| up 
to the rather general modelling of m ultip hase transfor- 
mation in a multi-component system 1171. All of these 
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models assume local equilibrium in the system, being 
consistent with the basic hypothesis of the classic irre- 
versible thermodynamics (CIT) 0, 0] . This assump- 
tion leads to the examination of a number of transport 
processes with small and moderate deviations from ther- 
modynamic equilibrium and, as a consequence, relatively 
slow movement of the interface can be predicted. In 
principle, such an approach can be extended to the case 
when the condition of equilibrium is violated locally at 
the interface, e.g., as it has been made for solute trap- 
ping and kinetic effects [2(j. However, the local equi- 
librium is missing both at the interface and within the 
bulk phases for rapid transformations such as rapid so- 
lidification 0| . In this case, the description of rapid 
phase trasformations might be provided by the formal- 
ism of extended irreversible thermodynamics (EIT) |22j 
which gives a causal description of transport processes 
and abandons the assumption of local equilibrium. An 
extension of the phase-field methodology for rapid trans- 
formation, which is caused by significant deviations from 
thermodynamic equilibrium, has been made recently |23| . 

The main purpose of the present paper is to describe 
a thermodynamically consistent model for rapid phase 
transformation in a binary system under local nonequi- 
librium conditions. Using the phase-field methodol- 
ogy, we derive the governing equations compatible with 
the macroscopic formalism of EIT and the microscopic 
fluctuation-dissipation theorem. 

The paper is organized as follows. In Sec. [HJ a ther- 
modynamic description of a considered system is given. 
We introduce the dissipative diffusion fluxes for heat 
and mass transport together with the phase-field rate of 
change as independent variables. In Sec. IIIII the gen- 
eralized Gibbs equation and entropy balance applicable 
to rapid advancing of the diffuse interfaces are given. As 
a starting point of the present phase-field model, an en- 
tropy functional is used in Sec. IIVI The analysis of the 
present phase-field model leads to the governing equa- 
tions for the hyperbolic system with dissipation. In Sec. 
IV1 a generalization of the hyperbolic phase-field model 
is given using the flux relaxation functions as well as a 
variational principle. In Sec. IVII the model equations 
are compared with the outcomes of the existing sharp- 
interface and diffuse-interface models. Finally, in Sec. 
I VI 1 1 we present a summary of our conclusions. 



II. DESCRIPTION OF THE SYSTEM 

A. Thermodynamic variables 

Let us consider an isobaric binary system at nonuni- 
form temperature T with no convective flow and given 
concentration of atoms A and B. The local equilibrium 
hypothesis establishes that the local and instantaneous 
correlations among the properties of the system are the 
same as for the whole system at a global equilibrium. 
Describing the nonequilibrium system as an ensemble of 



small local volumes in internal equilibrium, CIT |l9j is 
applicable to processes not too far from equilibrium. In 
addition to CIT, a local nonequilibrium formalism appli- 
cable to the strongly nonequilibrium systems has been 
developed in past two decades [22, 0, HH EH- As a 
phenomenological theory, this formalism is well-known 
as EIT [53, H3 which goes beyond the hypothesis of lo- 
cal equilibrium and avoids the paradox of propagation of 
disturbances with an infinite speed. 

A fundamental problem in attempting to describe sys- 
tems out of equilibrium is to select the relevant variables 
needed for a valid description of a nonequilibrium state. 
This problem has been intensively discussed in the liter- 
ature (see references in bibliographic overview A 
selection of the basic state space with the inclusion of 
the dissipative fluxes is formulated in EIT [24j and tested 
against experimental data |29| . Accordingly, we extend 
the classic set of independent thermodynamic variables 
by the inclusion of dissipative fluxes as additional basic 
variables. 

CIT is based on the local equilibrium hypothesis 0, 
[Tg| which assumes instant relaxation of fluxes to their 
steady-state values and describes the ensemble of atoms 
within local volumes by the Gibbs-Boltzmann statistics. 
In the standard formalism of the diffuse-interface using 
CIT, the set {C} of independent variables is assumed to 
consist of the conserved variables such as energy density 
e(r,t) and concentration X(r, t) — Xb/(Xa + Xb) of 
the B component in the system, and the non-conserved 
phase-field <I>(r, t) variable [where t is the time, and r is 
the position- vector of a point within system]. This can 
be expressed formally as follows: {C} = {e, X, $}. 

The extended space of independent variables E is 
formed by the union of the classical set {C} and the 
additional space {F} of the fluxes of heat q and solute 
J, and also the rate of change 9$ / dt of the phase-field 
variable, i.e. {F} = {q, J, d^/dt}. This yields 

E = {C} U {F} = {e,X,$} U {q,J,d$/dt}. (1) 

Here {F} is the space of fast non-conserved thermody- 
namic variables. 

There are, in fact, different possible choices of variables 
(fluxes in EIT, microstructural details in theories with 
internal variables), and the specific choice to be adopted 
depends on the aims of the description and on the prob- 
lems to be analyzed. This does not mean that different 
choices of variables are incompatible with each other. For 
instance, in the study of flowing polymer solutions one 
may select as independent variables either the viscous 
pressure tensor or the conformation tensor describing the 
average microstructure of the macromolecules of the sys- 
tem; a Legendre transform exists which allows one to 
pass from one description to the other, in a similar way 
as it is possible, in equilibrium thermodynamics, to pass 
from a description using internal energy as independent 
variable to a description using absolute temperature as 
an independent variable |30|. 
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Thus, our choice of the fluxes as variables does not ex- 
clude other possibilities. To justify our choice, we com- 
ment on the qualitative grounds, the meaning and the 
relevance of q, J, and d^/dt as variables. The fluxes 
q and J describe the exchanges of heat and matter be- 
tween the interface and the neighbouring bulk phases. 
The fluxes do not instantaneously follow by the classical 
Fourier and Fick laws, relating them with temperature 
and concentration gradients, it takes them some time 
(usually rather short) to reach the value predicted by 
the classical transport equations. Obviously, when the 
interface motion is fast enough, the delay effects in the 
dynamics of the fluxes may play a determining role. This 
happens, for instance, when the velocity V of the inter- 
face becomes comparable or higher than 1/t, I being the 
mean-free-path of the particles and r the relaxation time 
of the fluxes. Thus, in these circumstances, q and J be- 
have as independent variables with their own dynamics, 
which has important consequences on the dynamics and 
stability of the interface [31], |33 • 

The introduction of d$/dt as a further independent 
variable is motivated by similar, though slightly different 
consideration. Indeed, the space variation of $ is related, 
among other factors, to the width of the interface. Thus, 
including d^/dt as an independent variable allows for 
a more detailed description of the internal kinetics and 
shape of the interface. In the same way as in Newto- 
nian mechanics, where the initial position and velocity 
of a particle must be specified to determine their evo- 
lution, here we take both «3> and d&/dt as independent 
variables. If inertial effects are sufficiently low in com- 
parison with dissipative effects, d^/dt will be directly 
determined in terms of $ and its gradient by a dynamical 
equation. Otherwise, <!> and d$/dt will be independent 
and an equation for 9 2 $/9t 2 must be found. 

Consequently, with taking the above choice of vari- 
ables, one may distinguish between the two sets of inde- 
pendent variables as it follows. The variables from the 
set {C} are characterized as the slow variables, as their 
behavior is governed by conservation laws for energy and 
solute concentration plus the evolution of the phase-field, 
and as they decay slowly in time. In contrast, the inde- 
pendent space {F} consists of non-conserved variables 
with relatively high rate of decay. The variables from 
{F} differ from their classical value during intervals of 
the order of magnitude of the characteristic times t; for 
relaxation of the heat flux, solute diffusion flux, and rate 
of change of the phase-field variable. For time intervals 
much longer than these relaxation times Tj, the rate of 
variation of the fluxes can be ignored. 



solute can be considered as the times needed for smooth- 
ing of inhomogeneities of temperature and concentration, 
respectively, by diffusion. The time r$ of relaxation for 
the phase-field can be evaluated from the velocity of the 
diffuse interface moving through the local volume with 
the characteristic spatial length. Consequently, the rate 
of decay of the heat flux q, solute diffusive flux J, and 
phase-field rate of change d&/dt are estimated by the 
following characteristic times 



t t = a/V$, t d = D/Vi 



D- 



r* = l/V, (2) 



where a is the thermal diffusivity, Vt the finite speed for 
heat diffusion (i.e. the speed of propagation of temper- 
ature disturbances), D the solute diffusion constant, Vd 
the finite speed for diffusion (i.e. the speed of propaga- 
tion of concentration disturbances), V the velocity of the 
diffuse interface, and I the spatial length. 

For instance, the time tt is defined by the phonon- 
electron and phonon-phonon interactions for heat diffu- 
sion in metallic systems and it is estimated in Ref. 33] to 
be in the range of lCP 13 s < tt < 10~ n s. The time td is 
defined by the time for diffusion jumps of particles, which 
varies within a wide interval: I0~ 11 s < in < 10 s in a 
binary alloy system or inorganic solution [34j . In addition 
to this, the time r$ might be evaluated numerically from 
Eq. © assuming that the length I = Wq the width of 
the diffuse interface and the velocity V is the character- 
istic velocity for rapid adiabatic transformations. Thus, 
for numeric evaluation of r$ in a pure system, one may 
accept the following expression 



(3) 



where Q is the heat of transformation, \ is the heat ca- 
pacity (so that relation Q/x is considered as the charac- 
teristic temperature for adiabatic transformation), and 
Mo the coefficient for atomic kinetics. Taking the val- 
ues for pure nickel, e.g., Q/x — 418 K [35J, Mo = 0.52 
m/(s-K) [3(|, and Wq = 5 • 10~ 9 m, one gets r$ = 
2.30 • 10~ n s. This value for t$ fits well to the time of 
diffuse-interface kinetics which might be calculated from 
the "thin-interface" analyses of Karma and Rappel [33 
extended by Bragard et al. 13 . 

It is also reasonable to evaluate the relaxation time for 
the phase-field in a binary system using outcomes of the 
phase-field model via "thin-interface" analyses of Karma 
and Rappel (33 . Namely, for the nonisothermal solidi- 
fication of a binary system, Ramirez et al. |39l | derived 
the time t$ for the phase-field as a function of X and 
It is described by 



B. Relaxation times 



7$ = 



r 



1 , Wo 

h a\ai- 

Mo 



D 



DQ m(l-k)X 
ax + 1 + k - (1 - fc)<& 



■(4) 



Generally, the relaxation times Tj represent physically 
reasonable time estimations for the spontaneous return of 
the system to the steady state after some sudden pertur- 
bation. The relaxation times tt and td for the heat and 



For numeric evaluation, we accept the following ma- 
terial parameters for a Cu-Ni metallic system in Eq. 
(@J: diffuse-interface width W = 1 • 10~ 9 m, Gibbs- 
Thomson coefficient T — 1.3 • 10~ 7 K-m |4jj, atomic 
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TABLE I: Relaxation time for the fluxes of heat, solute diffusion, and phase-field 



System tt (s) td (s) t<s> (s) 

Carbon tetrachloride 2.50 • 10" 13 , Ref. [42] 

Benzene 1.22 ■ 1(T 13 , Ref. [42] 

Nickel 1.20 ■ lO -11 , Ref. [42] 2.30 ■ 10 -11 , Ref. [44] 

Diluted alloy Ni - 0.7 at.% B 1.54 ■ 10" 11 , Ref. [45] 

Concentrated alloy Cu - 30 at.% Ni 0.75 ■ 10~ n , Ref. [46] 



kinetics coefficient /io = 0.24 m/(s-K) 40], constants 
ai = 0.8839... and a 2 = 0.6267... |37|, solute diffusion 
constant D = 3 • 10~ 9 m 2 /s [4jj, thermal diffusivity 
a = 1.5 • 10~ 5 m 2 /s adiabatic temperature (rela- 
tion of latent heat and heat capacity) Q/x — 402 K [4l| . 
slope of the liquidus line m = 4.38 K/at.% solute 
partitioning coefficient k = 0.81 |4l|. As a result follow- 
ing from Eq. (0J, one gets r$ = 7.92 • 10~ n s for the 
values of X = 70 at.% and $ = 0.5. 

The values for the relaxation times for some pure and 
binary systems are summarized in Table [I] It can be 
seen, e.g. for metals and alloys, that even though the 
heat speed Vr is much larger than the solute diffusion 
speed Vd , the relaxation times for q and J have the same 
order of magnitude, i.e. tt ~ To. Therefore, a front of 
the heat profile moves with a speed much higher than a 
front of the solute diffusive profile. However, due to the 
fast thermal diffusion, a >> D, the relaxation of the heat 
flux q proceeds approximately at the same characteristic 
time as the relaxation for solute diffusion flux J. 



III. ENTROPY APPROACH 

A. Generalized Gibbs equation 

For the local nonequilibrium system described in Sec. 
Inl we postulate the existence of a local generalized en- 
tropy density s whose set of variables is the extended 
space E by Eq. QJ. The generalized Gibbs equation for 
s is described by 

ds(e, X, $, q, J, d<S>/dt) = ds e (e, X, $) 



+ds ne (q, J, d*/dt) = ^de + ^dX + 

oe dX o<P 



clq H — • aJ ■ 



ds 



dq 



dJ 



d 



d{d$/dt) V dt 



(5) 



In Eq. J5J), s e is a local equilibrium contribution defined 
on the set {C} of the classic slow variables {e, X, $}, and 
s ne is a flux-dependent purely nonequilibrium part of the 
generalized entropy defined on the space {F} consisting 
of the fluxes {q, J, d$>/dt} as the independent fast vari- 
ables. 



The derivatives of the entropy density with respect to 
classical variables and their fluxes appearing in Eq. (JSj 
are described by 

ds e 1 ds e A/i 
~de = T' ~dX = T~' 



ds e _ v ,ds A 

9$ - {1 ~ X} ^¥ 



X 



ds B 
d<$> ' 



ds ne 
dq 



ds 



~a q q, 



dJ 



-a, J, 



ds r . 



9$ 
dt' 



d(d$/dt) dt ' (6) 

where A/i = [la — /is is the difference of the chemical 
potentials for components A and B, respectively, and 
sa and sb are the entropies for pure components A and 
B, respectively. The chemical potentials and entropies 
of components can be chosen for every concrete system 
(see, e.g., Refs. 

In Eqs. JBJ, the coefficients ai are scalars which do not 
depend on q, J, and 9$ / dt and are assumed to be 



Lr 2 ) 



cv, 



TD \ dX J 



Q!0 



a - 



t*WqQ 



(7) 



T,X 



where k is the thermal conductivity, ao a dimension- 
less factor (dependent on the model of the diffuse in- 
terface, specifically leading to the sharp-interface asymp- 
totic limit), and Q the heat of the transformation. 
After integration, Eq. (JSJ) can be written in the form 

s(e, X, $, q, J, d^/dt) = s e (e, X, $) + s ne (q, J, d$/dt), 



s ne {q,J,d<S>/dt) 



a, i 



-q ■ q 



-J- J 



a<t, 
~~2 



9$ 
~dt 



(8) 

Consequently, we arrive to a generalized entropy density 
given by an expansion around its local equilibrium value 
up to second-order in the fluxes. In the limit of infinite 
speeds (Vr — * oo, Vd —* oo, and V — » oo), one gets 
tt — > 0, td — > 0, and r$ — > 0. In such a case, the 
term s ne vanishes and Eq. (JSJ gives the entropy density 
s e (e,X, $) for local equilibrium system. 
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B. Entropy balance 

For the system described by the extended set E of vari- 
ables, Eq. QJ, the local balance laws for the energy and 
concentration are given by 



* + v. f -o, § + v./=o, 

and the evolution of entropy density is defined by 
ds 



at 



V • J s = crs- 



(9) 



(10) 



The change of the total entropy S in time t is described 
by 



dS_ _ fdS_ 
dt ~ \dt 



dS 
~dt 



(11) 



where 



\dl) = ~ J V X 7 ■ J ' sdv = - J Js-ndw, (12) 

is the external exchange of entropy due to entropy flux 
Js and 



dS_ 

dt 



asdv, 



(13) 



is the internal production of entropy due to dissipation 
within the system. In Eqs. (|12|) and l|13fl : co is the outer 
surface of sub- volume v, n the normal vector to the sur- 
face, and as the local entropy production. 



plicit evolution equations for the variables including the 
relaxation terms. 



A. An entropy functional 



Now we use an entropy functional of the following form 



S= j [s(e,X,$,q,J,d$/dt) 

J V 



-^\Ve\ 2 -f\VX\ 2 -f\V<i>\ 2 ]dv. (14) 



Here e e , e x , and are constants for the energy, concen- 
tration, and phase-field, respectively. In the functional 
the gradient terms |Ve| 2 , VX\ 2 , and |V$| 2 are used 
to describe spatial inhomogeneity within the fields ac- 
cording to previous diffuse-interface models [a, lift 
It is logical to include gradient terms in Eq. (|14fl [of the 
so-called "Ginzburg-Landau form"] because, as stressed 
before, our interest is focused on interfaces with steep 
gradients. In addition, the extension Q gives the en- 
tropy density s based also on the fluxes q, J, and d$/dt 
as independent variables. 



IV. HYPERBOLIC PHASE-FIELD MODEL 

In this section, the important class of hyperbolic mod- 
els with dissipation is considered. We work out the ex- 



To obtain the evolution of the entropy, Eq. l(TT|l . 
and consider the several parts of the entropy exchange, 
Eqs. l(T2* )l -l(T3 )l . we differentiate Eq. (JT3J with respect to 
time. Combining the terms, after some algebra one ob- 
tains 



J 



dS_ 

~dt 



ds 



elV 2 e 



dv 



ds 
dX 



elV 2 X 



dX 
~dt 



dv 



ds 
9$ 



+ 



ds ( dq\ ds I dJ 

d~q\di) + a/i W . 



dt 



d{d$/dt) \ dt 2 



dv 



de 



dt 



dX 



< V„e + si — V„X + ei — V„$ 



dt 



d<P 



dt 



du, 



(15) 



where V„ is the gradient vector pointed by the normal vector n. 
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Now we substitute the balance laws for energy and 
concentration, Eqs. into Eq. I|15|). and then use the 



theorem of divergence. One gets 



dS 
~dt 



V„e - 



ds 
de 



+ e 2 e V 2 e q 



dX 
~dt 



Vr,X 



ds 
dX 



si\7 2 X J, 



d<f> 



V„$ >duj 



g-V 



ds dq -• „ 



9s 



5s <9J a$ 



2l + e 2 v 2 $ 



<9 2 $ 



d(d$/dt) dt 2 



dv, 
(16) 



where g„ and </„ are the diffusion fluxes pointed by the 
normal vector n. 

Using Eq. J§}, the change of the entropy Eqs. 1111) - 
(|T3j) . is obtained from Eq. (JTSJ. This yields 



a s dv, 



where 



Js=e 2 e 



V„e - 



ds 
d~e 



e 2 e V 2 e ) q n 



* Vdt 



(17) 



(18) 



is the projection of the entropy flux vector on the normal 
vector n, and 



a s = q ■ 



B. Governing equations and thermodynamic 
consistency 

Relation (|18f) is well known from the phase-field model 
based on CIT (see, e.g., Ref. whereas the entropy 

production (|T1?)l includes the additional terms —a q dq/dt, 
—OijdJ/dt, and —a^^/dt 2 related to the nonequilib- 
rium part of the generalized entropy. This is due to the 
special form for entropy Eq. (JSJ , and has a clear physical 
meaning: far from equilibrium, the dissipative fluxes pro- 
vide ordering that leads to a decrease of the entropy pro- 
duction near a steady state as compared with the local- 
equilibrium state characterized by the same values of e, 
X, and $. 

The production as of the generalized entropy, Eq. © 
is positive due to the statement of the second law of ther- 
modynamics. This condition implies a relation between 
fluxes and forces which, in the simplest cases, is assumed 
to be linear. For Eq. I|19|) . it gives the following set of 
equations: 

- evolution equations for heat and solute 



usion fluxes 



(1 
J 



(M){ 



ds 
~de ' 

ds 



r 2 V 2 e 



dq 
'dt 



v[ dx + ^ x 



a, 



dj_ 

~dt ) 



evolution equation for the phase-field 



(20) 



+J- 



_ , ds ~ , \ d J 



where 



<9 2 $ 



<9$ 

"dt 



ds - - d 2 & 



> o 



(19) 



is the local entropy production which has a bilinear form 
in the fluxes (q, J, and d&/dt) and their respective con- 
jugate forces (expressions inside the square brackets). 



(M) 



( M KR M R , 



M™ M r . 



(22) 



is the matrix of mobilities for thermal and solutal trans- 
port, and is the mobility of the diffuse interface. 
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Dependent on composition, the interface mobility is as- 
sumed to be 



Ma, 



(l-X)Mt + XMZ >o, 



(23) 



where M$ and are the interface mobility for the 
transformation in pure systems consisting of A or B 
components, respectively. In various formulations of the 
phase-field model [2(J|33, the mobilities of Nl£ and 
are proportional to the atomic interface kinetic coefficient 
Ho and inversely proportional to the interface width Wo, 
so that Ma, ~ Hq/Wq. 

The matrix lll'L'1) of transport and the interface mobility 
(|23J) are assumed to be positively defined for the positive 
entropy production as- The matrix can be consid- 
ered as symmetric, so that the matrix can be regarded as 
being positive with the inequality: M ee M xx > M 2 X . Note 
that the linear phcnomenological laws given by Eqs. 12U|) 
and (|21|l adopt the representation theorem of isotropic 
tensors |48) according to which fluxes and forces of dif- 
ferent tensorial rank do not couple as far as linear re- 
lations arc involved (this independence of processes of 
different tensorial rank is also known as the Curie princi- 
ple) . In our case, the vectors of heat and solute diffusion 
fluxes cannot give rise to the flux of the scalar phase-field 
flux in a linear description. More complicated nonlinear 
relations between fluxes and forces consistent with posi- 
tive entropy production in EIT are considered elsewhere 

For simplicity, we ignore both kinds of "cross coupling" 
effects in Eq. flU, so that M ex = M xe = 0. Then, 
substitution of the fluxes from Eq. H2U|) into the balances 
©, respectively, gives 

- the governing equation for energy density 



d^e 
dt 2 



de 
dt 



= -V • 



M fifi V 



ds 
de 



+ elV 2 e 



the governing equation for solute concentration 



d 2 X dX 
di? + ~dt 



= -V- 



M XX V 



( ds 
(dX 



el\7 2 X 



(24) 



(25) 



in which tt = a q M ee is the relaxation time for the heat 
diffusion flux, and td = ajM xx is the relaxation time for 
solute diffusion (see Eqs. J5J) and Table HJ). After simpli- 
fying the transformation, Eq. H21|) leads to 
- the governing equation for the phase-field 



9 2 $ 9$ ^ (ds 2r7 _ 



(26) 



where r$ = a^M^ is the timescale of the phase-field ki- 
netics. According to Eq. (|26|l . the acceleration d 2 &/dt 2 
of the phase-field appears due to introduction of both 
<& and d^/dt as independent variables and characterizes 
inertial effects inside the width of diffuse interface. 

Equations l|24|) - (|26|l are the central outcome of our pro- 
posal [or, to mention a more complicated setting, we 



could also refer to equations (|20|I - H21() ]. The role of the 
relaxation times is clear: they characterize the delay with 
which q and J reduce to their classical expressions (corre- 
sponding to classical transport equations) , and the delay 
with which the inertial effects in the dynamics of the 
interfacial region are lost. The relaxation terms may 
be neglected in many circumstances, but become crucial 
in some important situations, leading, for instance, to a 
maximum possible value for the speed of advancement of 
the interface (in contrast to classic theory which allows 
for an infinite speed of propagation) , and to the possibil- 
ity of oscillatory phenomena in the width of the interface. 
Thus, the role of the new terms is not simply to add some 
new undetermined parameters (relaxation times) allow- 
ing for an improved fit of experimental results. These 
terms also play an important conceptual role, as they 
drastically change the possible kinds of behavior of the 
system. 

Some comments on the consistency of our proposal can 
be outlined. First of all, we may refer to its internal con- 
sistency as a thermodynamic (macroscopic) theory. Sec- 
ond, one must check its consistency with microscopic de- 
scriptions based, for instance, on kinetic theory, or on lin- 
ear response theory, or in other statistical (microscopic) 
theories. Finally, one must check its consistency with 
experimental results. 

Here, we comment on the internal thermodynamic 
consistency and, in the next Section, we shall refer to 
its consistency with a statistical theory, based on the 
fluctuation-dissipation theorem. In this theoretical pa- 
per we do not refer to experimental results. We assume 
that a consistent noncquilibrium thermodynamic theory 
should satisfy two main conditions: 

(i) the generalized or extended entropy must be maxi- 
mum at equilibrium; 

(ii) the entropy production must be positive. 

To these two conditions one could add two more re- 
quirements: 

(iii) the second differential of the entropy with respect 
to its basic variables (which is related to the dynamics 
of the variables) must be negative in order to lead to dy- 
namically stable solutions; 

(iv) the generalized equations of state obtained by differ- 
entiation of the generalized entropy must have a physical 
meaning consistent with experiments. 

It can be seen immediately that the essential condi- 
tions (i) and (ii) are satisfied in our proposal. Indeed, 
the form © and (|I4|I of the entropy guarantees that ho- 
mogeneous equilibrium state has the maximum entropy 
as compared to nonequilibrium states with the same local 
values of e, X and Furthermore, introduction of the 
constitutive equations H2U|) - (|2I|I into the expression i|19|) 
of the entropy production yields for the latter a definite 
positive expression: 



as = (V, j) (My 1 + M7/ 1 > 0. (27) 
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As we noted, the transport matrix (A4) and the inter- 
face mobility are assumed to be positively defined 
for the positive entropy production, as > 0. If one 
included higher-order nonlinear terms into the entropy 
(JSJ or in the constitutive equations l|20|l - H21() . thermody- 
namic consistency would be more difficult to check than 
in our second-order approximation (jSJ). This approxima- 
tion is sufficient to deal with a wide range of physical 
problems. 

We shall not deal with conditions (Hi) and (iv), which 
are subtler and typically involve nonlinear effects. For an 
indication of their analysis in some situations involving 
only q as nonequilibrium variables, the reader is referred 
to the monograph 24] . 



V. GENERALIZATION OF THE MODEL 

The governing equations (|24l) - (|26|l present causal prop- 
agation of heat and mass signals and a dissipative-wave 
advancing of a diffuse interface. We generalize them into 
evolution equations which are nonlinear in time. First, 
equations of state are considered from the point of view 
of the relaxation functions for the fluxes. Second, non- 
linear evolution equations of a general type are derived 
from a variational formulation. 



A. Relaxation functions for the fluxes 

Let's take into consideration a prehistory of the change 
of the phase-field in a point of a system. Such a prehis- 
tory must be taken if the system is not in local equi- 
librium. We shall use a functional description with a 
memory function. 

We use the entropy functional (|T4*jl . as before, to derive 
the equations of the model. In the absence of local equi- 
librium, one may incorporate the prehistory of the diffu- 
sion process. Then, the connections between the fluxes, 
q, J, and d^/dt, from the one side, and driving forces, 
V(5S/Se), V(5S/SX), and 5S/5&, from the other side, 
are defined by the following integral forms: 
- relaxation of the heat flux 



q(r, t) 



I D q (t-t*)W 

J — GO 



SS{t*,r) 
6e 



dt* 



(28) 



relaxation of the solute diffusion flux 



J(r,t) 



D 3 (t~t*)V 



SS{t*,r) 
SX 



dt*, (29) 



relaxation of the phase-field rate of change 

1 9mt) ~ ^ D^t-n^A^ (30) 



M dt 



6$ 



where Dr = {D q , Dj,D^} are the relaxational kernels for 
the fluxes, and the variational derivatives are obtained 



from the following expressions 

SS ds 9 o SS ds 
— = 1- e V e = 

Se de e ' SX dX 



elV 2 X, 



S_S_ 

~sl> 



ds 



(31) 



After substitution of expressions for the heat flux re- 
laxation, Eq. 1)28(1 . and the solute diffusion relaxation, 
Eq. (|29|l , into the balance laws for energy and solute con- 
centration, Eq. (JSJ, respectively, one can get the following 
integro-differential equations 



de(r,t) 
dt ' 

8X(r,t) 
dt 



Se 



= -V • 



D 3 {t-t*)V 



SS(t*,f) 
SX 



eft* (32) 



Together with relaxation of the phase-field, Eq. H30|) . the 
general system evolution during phase transformation is 
described by Eqs. (|3^|) . 

When the relaxation functions Dr are specially de- 
fined, Eqs. ptljl and l|3"2*|l can be reduced to known mod- 
els. Particularly, for the important class of dissipative 
and hyperbolic models, one can take the relaxation ker- 
nels in the following forms 



Dr(0) = const, 
D R (0)5(t-t*), 



D R (0) exp - 



t-V 



wave propagation, 
dissipation, 

wave and dissipation, 
(33) 



where D R (0) = {D q (0), Dj(0), D${0)} are the relax- 
ational kernels for the fluxes at present time t — t* , and 
t = {tt,td,t$>} are the characteristic relaxation times 
for the fluxes. 

The different transformations within the diffuse- 
interface are described by different kernels in the inte- 
grals (I2H1>- (EOI - As it follows from Eq. (J33J), the relax- 
ation functions Dr describe the memory of the system 
by assigning different weights to different moments in the 
past. Dissipation corresponds to a zero-memory transfor- 
mation, i.e. the only relevant contributions are the "last" 
ones. In contrast to this situation, the infinite mem- 
ory transformation with Dr = const leads to undamped 
wave propagation of the heat, solute, or the interface 
advancement. In between, the combination of the wave 
and dissipative regimes described by the exponentional 
law can be observed during rapid phase transformations. 
This is the case of hyperbolic phase- field model described 
in Sec. IIVI For the latter case, the relevance of all con- 
tributions to the fluxes decreases as the system moves to 
the past. 

In Sec. IIVI the model macroscopic consistency of the 
statements of EIT has been shown. Now, the consistency 
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of our macroscopic approach with a microscopic descrip- 
tion is verified in relation to the outcomes following from 
the fluctuation-dissipation theorem. 

The memory functions introduced in Eqs. l|28[l - i|3U|) 
may be related to our analysis of the dynamics of the 
fluxes q and J and of d<fr/dt proposed by constitutive 
equations l|2U|) - (|21|l . To do this, first, we may consider the 
fluctuation-dissipation theorem relating response mem- 
ory functions to the time-correlation function of the cor- 
responding fluxes (see, e.g., Ref. ^3). This will allow us 
to show the consistency of our macroscopic formulation 
with the microscopic basis provided by the fluctuation- 
dissipation theorem. 

The corresponding expressions are 



which may be rewritten as 



D q (t-t*) 
Dj(t-t*) 



1 



k B T 2 
1 



WW**))** 



knT 



(J(t)J(t*)) eq , 



D*(t-t*) = —-(d t $(t)dMt*)) eq . 



(34) 



Here k B is Boltzmann's constant, q, J, and 9t$ stand 
for the microscopic operators for the heat flux, diffusion 
flux and the time derivative of <E>, respectively, and (...)eg 
means an average over an equilibrium ensemble in statis- 
tical mechanics (as, for instance, the canonical one). 

Relations (|34|) play an important role in modern statis- 
tical mechanics, and may be formally derived from the Li- 
ouville equation in the framework of linear-response the- 
ory or from information theory |49ll50|. However, from a 
practical point of view, the computation of the evolution 
of the microscopic operators for q, J or <9 t $ on purely mi- 
croscopic grounds is an overwhelming task exceeding ac- 
tual capabilities. Such an evolution is either obtained by 
computer simulations, or tentatively given by a reason- 
able form inspired on phcnomcnological grounds. Thus, 
our evolution equations (|2UI) - (|21[I for q, J and 9t$ may be 
considered as a macroscopic modelling of the evolution 
of the fluxes, which according to Eq. (|34|l is equivalent 
to proposing a form for the corresponding memory func- 
tions introduced in Eqs. lf^S|) - (|3l7|) . In general terms, it 
could be said that, according to Eq. ffi^i . the study of 
the evolution of the fluxes around equilibrium is equiva- 
lent to the determination of the corresponding memory 
functions. 

Constitutive equations 120() - (|2I() imply that fluctua- 
tions of q and J near a homogeneous equilibrium state 
will decay exponentially as q(t) = q(0) exp(— t/rx) and 
J(t) — J(0) exp(— t/T£>). Introducing these expressions 
into Eq. J3H one obtains 



D q (t-t*) 



1 ^ - ( t -t* 

'<f(0)<f(0)) eg exp 



k B T 2 
I 



— <?(0)?(0)) eg exp \-~- ) i ■->'') 



D q (t~t*)=D q (0)exp 
D j(t -**) = £>j (0)exp ( - 



t - t* 

t - t* 
TD 



(36) 



Indeed, when the microscopic expressions for q and 
J corresponding to an ideal gas are introduced into 
Eq. (|36[) and the equilibrium average is performed (over 
a Maxwell-Boltzmann distribution function), the results 
for the D q (0) and Dj(0) are equivalent to those obtained 
from the kinetic theory of gases in the time-relaxation 
approximation |49l |. 

Note, finally, that the usual transport coefficients 
(thermal conductivity, diffusion coefficient) may be ob- 
tained (when the relaxation time is sufficiently short) by 
integration of Eq. ifMj) . as 



A 



D 



1 



k B T 2 J_, 
I 



knT 



(q(t)q(0))e q dt, 



(J(t)J(0)) eq dt, 



(37) 



which are the well-known Green-Kubo formulae for trans- 
port coefficients 0, 0, Thus, our macroscopic 
formalism is consistent with the microscopic fluctuation- 
dissipation theorem. It provides, in fact, a phenomeno- 
logical complement to the fluctuation-dissipation expres- 
sions, which are the formal expressions from which it is 
difficult to obtain on exact grounds the form of the mem- 
ory functions. 



B. A variational principle and Euler-Lagrange 
equations 

We assume, as above, that the generalized entropy 
density s is a continuous and differentiable function de- 
fined by the local equilibrium contribution s e and flux- 
dependent nonequilibrium part s ne with the total set of 
variables JIJ and generalized Gibbs equation JSJl. The 
balance equations for the heat and solute are the same, 
Eqs. ©, and the local evolution of the entropy density 
is described by Eqs. Q1U[I. 

A generalization can be given by introducing the gen- 
eralized terms for derivatives into the entropy density 
with respect to classical variables (e, X, $) and fluxes 
(q, J, 9$ / dt) , and also by introducing general forms of 
the entropy flux J$ and the source as in Eq. i|IU|) . De- 
pending on their own tensorial character, these are 



ds 

ds 
9$ 



ds 
dX 



— 



— =f3 e 2 (e,I q )q, 
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0s 

dJ/x 
( — \ 



0$ 

01^,1^)—, 



J s = Pi(e, I q )q + ff{X, Ij)J + %)^, 



a s =l3t(e,I q ) + 0? (X, I, ) + (3f (<&, % ) , 



where 



I q = Q- 9, 



I j = J-J, 



Io<s> 
ST 



d<S> 
0t 



(38) 



(39) 



are the single scalar invariants of the extended set (JTJ 
of variables, and j3i are the scalar functions depending 
on classic variables (e,X, $) and invariants Ii. Then, 
utilizing Eqs. I|38|) . the generalized Gibbs equation (J5J 
gives the time derivative of the entropy density as follows: 

! -«(.,/.)!+«(., w-f? 

+ft V (A',/ J )^+ft v (X,/ J )J-^ 

9$ 9$ 2 $ 

+/f(*.¥ft+^»%%F' (40) 

Locally, Eq. (|lt)|l is satisfied as a balance law and, for 
the entire system, one can postulate extremal condition 
in the Lagrangian form of£ = f v (ds / dt+V ■ Js—&s)dv — > 
extr, implying an extremal difference between the "ki- 
netic" part / (ds/dt + V ■ Js)dv and the "potential" part 



J asdv for the whole nonequilibrium system. Then, the 
entropy density satisfies the following variational princi- 
ple p| 



SC = :) I ,lr [ |i 



V • J s - <J S }= 0, 



(41) 



in which the variation S is carried out only on the noncon- 
served flux variables q, J, and d$/dt, i.e. S is taken only 
over the space {F} from the set (QJ while the variables e, 
X, and $ from the set {C} remain constant during the 
variation. Also, during the variation, the tangent ther- 
modynamic space [time and spatial derivatives from the 
set iJIJ] is fixed. From this it follows that Eq. (|41|> is a 
variational principle of a restricted type. 

Using balance laws Q, substitution of Eqs. 1)38(1 and 
(|40|1 into variational principle (|41|l leads to 



dv 



dJ 



J* 



J - flf 



A 0$ * 0$ / *0 2 $ *\ 9$ 



= 0. 

(42) 



Variation of Eq. ((421) is obtained by taking as constants 
the time derivatives, gradients and divergences. Since 
SI q = 2q-Sq, SL = 2J-SJ, and Wo* = 2(0$/0tW0$/0t) 

____ dt 

from Eq. 142(1 one gets 



dv 



3 « W«v.» + «« 



81a dl, 



dt 



Bin 



dt 



V/3 3 e 



2 f* 

Ol a 



8q 



dv 



dv 



x 



dlj dl 



dla® 
ST 



0$ 
at 



J(V- J)+/3& + 2 



0j 



.0/3 



0J 



St 



0/?f SJE Sg 



s/ 6 



ST 



sf JJ '^ + VA 



4 s^ 

^ 2 St 2 



J s/ ? 



2 0^s^/s$V 

dlo* dt 2 \dt f ' 3 dlo* dt 



ST 



TIF 



0$ 
0t 



= 0. 



(43) 



Due to arbitrary variation of Sq, SJ, and 5(0$/0t), the 
Euler-Lagrange equations directly follow from Eq. ((43(1 . 
These are 



- evolution equation for the heat flux 



( 0/3f e \ 0? 
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~{d(3i 


d[3f\ 




91 J 



V-q 



dJ a 



1 



V/3 3 e , (44) 



evolution equation for the solute diffusion flux 



dL 



J J 



df%_ 
dL 



x 



an 



dJ_ 
' ~dt 

V • J 



x 



m 

dL 



(45) 



evolution equation for the phase-field 



df3? 




Bio® 

v TTT 



3 • 



(46) 



where W is the unit tensor of second rank. 

Eqs. H44|) - (|46|l are nonlinear evolution equations for q, 
J, and d&/dt and they are of the general form of evo- 
lution equations l|2U|) - (|21|l . Indeed, the nonlinearity is 
clearly seen from the following form of these equations: 

T T (e, q)^+q = M ee (e, q)Vp e 3 , 
t d (X, J)— + J = M XX (X, J)V/3f , 

where and Mi are the functions of the classic set 
{C} — {e, X, $} as well as nonlinear functions of the 
fluxes which can be explicitly found from Eqs. i|44|l - i|4t)|) 
and relations i|38|) - (|39|) . Thus, taking the generalized evo- 
lution of the entropy density, Eq. i|40|) . and using varia- 
tional principle l|41|) . we arrive to nonlinear general evo- 
lution equations for fluxes, Eqs. (|47|l . which might be 
merely reduced to the evolution equations l|2U|) - (|21|) of 
the hyperbolic phase-field model. 



VI. RELATION TO EXISTING MODELS 

It is interesting to note that sharp-interface and 
diffuse-interface models with relaxation of fluxes have 
been used to describe transient processes in various 
nonequilibrium systems (see monograph [53, Chapter 4). 
Therefore we synthesize here several previous and very 
recent results in comparison with the developed hyper- 
bolic model (Sec. llVfl and generalized model (Sec. 01 of 
rapid phase transformation. 



A. Superconductivity 

Ginzburg and Landau established their variational 
principle for the continuous transition from the normal 
to the superconducting phase @. They used a free en- 
ergy density with a gradient term which has been further 
used in many phenomena (e.g., in description of spinodal 
decomposition [53j or crystal growth [5J|). As a logi- 
cal extension, the transition between the normal and the 
superconducting phases can be described with the de- 
lay given by equations of the hyperbolic model [starting 
from the functional of the form (|14l) ] or using generalized 
models with memory, Sec. IV Al 

Generally, equations l|24|) - (|26|l are consistent with the 
generalized entropy density given by Eq. (JSJ). The equa- 
tions are reduced to the classic equations from Refs. 
|sL l53l |54| when the times tt, td, and r$ tend to zero. 
Furthermore, the entropy density © together with the 
evolution equations l|20l) has been justified microscopi- 
cally [24L Efij for the one-component system and from 
Grad's procedure for monatomic gases. 

The choice of thermodynamic potential is important, 
as it governs the transition from metastable state to the 
stable one. Normally, the potential for transition is in- 
cluded in the expression for entropy density (or for free 
energy density) in the form of a double-well function 
or by a monotonically increasing function incorporating 
nonequilibrium conditions at the interface |8|, Il2l |38( . In 
the present paper, we do not give an explicit form of s e 
in Eq. © and present governing equations l|2"4l - (|27)|) [or 
variational derivatives (|31f) ] in a general form. The choice 
of the thermodynamic potential might be given for the 
problem under consideration. 



B. Glass transition, structural relaxation and 
phase separation 

Jackie et al. .55] considered isothermal phase trans- 
formation in the presence of additional slow structural 
relaxation variables. Considering the dynamics based on 
the relaxational chemical potential, these authors refer 
their model to systems with phase separation and to slow 
structural relaxation in polymeric solutions in the prox- 
imity of the glass transition temperature. The calcula- 
tion has shown that, even at the early stages of phase 
separation, equation for chemical potential with memory 
may give pronounced deviations from the predictions of 
classic linear Cahn-Hillard's model |H^ . 

Phase separation during spinodal decomposition may 
proceed under local nonequilibrium conditions in solute 
diffusion field offered by rapid quenching. As it has 
been demonstrated in computational modeling |56| . the 
rapidly quenched liquid mixtures under decomposition 
exhibit non-equilibrium patterns, evolving with univer- 
salities different from those extracted from the Cahn and 
Hillard's model. 

Local nonequilibrium separation in liquids can be de- 
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scribed in terms of EIT as a model for isothermal spin- 
odal decomposition in a binary system [23j , in conditions 
of large deviation from thermodynamic equilibrium. The 
dynamics of the diffusion flux J as a fast variable from 
the set JQ is consistent with the characteristic time of 
local rearrangement of particles (atoms or molecules) or 
with the time of relaxation of diffusive flux to its local 
equilibrium steady-state value. The model equation for 
spinodal decomposition of a binary system is the gener- 
alized Cahn-Hillard equation of the form of Eq. (|25[1 for 
local nonequilibrium solute redistribution. In this case, 
the dynamics of rapidly quenched decomposition is de- 
scribed for short periods of time or large gradients of 
chemical composition. 



C. Shear flow, viscoelastic fluids and 
diffusion-reaction systems 

The system of coupled evolution equations l|2(J|l - H21|) 
describes, in fact, a process of phase separation under 
shear if temperature is replaced by viscous pressure ten- 
sor. In this case, one may get the condition defining the 
spinodal line in non-equilibrium states [see monograph 
|30|. Chapter 6]. As a reduced one, equation of type l(2T]|) 
or l|32|) . endowed with homogeneous Dirichlet boundary 
conditions, has been introduced to model the behavior of 
certain viscoelastic fluids and to predict the velocity of 

flow H3. 

In addition, equation of type (|S2*|l is used to predict 
a wavefront in time-delayed reaction-diffusion systems of 
the generalized Fisher's equation [5^|. The speed of the 
travelling wave depends on the relaxation time and there- 
fore spreading of population in reaction-diffusion system 
might be predicted with great flexibility. One of the con- 
sequences of this equation, reduced to Eqs. (|2*1|) or lt2*5|) 
for modelling a hyperbolic reaction-diffusion system [with 
e e = or e x = 0, respectively], might be considered in 
an exciting example suggested by Fort and Mendez in 
Ref. |59|| for neolitic advancing of human groups across 
Europe. They have shown, in particular, that hyperbolic 
reaction-diffusion equations of type of Eqs. (|24H or 125|) 
predict the population spreading during the European 
past in agreement with the existing archaeological data. 

D. Rapid solidification 

At deep supercoolings in a solidifying system, or at 
high velocity of the solid-liquid interface, it is necessary 
to take into account local nonequilibrium effects in so- 
lute diffusion and to use a non-Fickian model for trans- 
port processes which is compatible with EIT pH |40| . 
The problem of rapid solidification within the sharp- 
interface limit is described by generalized Stefan prob- 
lem (so-called "self-consistent hyperbolic Stefan prob- 
lem" [i^], ) which takes into account local nonequilib- 
rium both at the interface and within the bulk phases. In 



such a case, the spatio-temporal evolution of solute con- 
centration is described by the partial differential equation 
(|25|1 of a hyperbolic type [with e x = 0] which takes into 
account the relaxation of solute diffusion flux to the lo- 
cal thermodynamical equilibrium in a rapidly solidifying 
system. 

Advancing of the diffuse-interface with a higher ve- 
locity comparable with the solute diffusion speed is also 
described by t he p hase-field model with relaxation of the 
diffusion flux [23]. It has been shown that choosing the 
concrete form of the entropy (as the thermodynamic po- 
tential), one may recover the existing models based on 
the CIT and analyze solidification under local nonequi- 
librium conditions. 



E. Motion of antiphase domains 

In the description of diffuse interface kinetics, Allen 
and Calm 10] proposed a model for evolution of the 
non-conserved order field during antiphase domain coars- 
ening. For isotropic interfaces, gradient flow gives the 
Allen-Cahn equation by taking r$ = in Eq. This 
equation is true in the case of low inertial effects in com- 
parison with the dissipative effects. With finite relax- 
ation time, r$, and finite acceleration, d 2 Q/dt 2 , Eq. I|26|) 
predicts the evolution of coarsening with relaxation. It is 
reasonable to say that the generalized Allen-Cahn equa- 
tion (|26|l is true for the case of significant inertial effects 
during the motion of antiphase domains. 

As an advancing of the Allen and Cahn's model, the 
process of the interface motion by mean curvature with 
delayed response has been analyzed recently. Rotstein 
et al. |(3(j developed the phase-field model based on 
equations similar to Eqs. (|28|l and i|3U[l . These authors 
described the first-order transition with the delayed re- 
sponse of the system due to slow relaxation of internal 
variables. Using the exponential relaxation function for 
wave and dissipative mode, Eq. which leads to the 
hyperbolic phase-field model, the dynamics of the per- 
turbed motion of interface by mean curvature has been 
considered. It has been shown in Ref. |60| that internal 
relaxing effects induce damped oscillations in the inter- 
facial motion during crystalline coarsening. As opposed 
to the classic parabolic phase-field model, the hyperbolic 
phase-field model predicts these interfacial oscillations in 
qualitative consistency with the oscillations on the sur- 
face of quantum crystals 61] and in crystallization waves 
in helium |62|. From a mathematical viewpoint, a search 
for existence and uniqueness of the solution and some 
well-posedness results for the problem of motion by mean 
curvature using the phase-field model with memory are 
beginning to be presented (see Ref. [d^I. 
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F. Complex (dusty) plasmas 

Recent investigation into the field of complex (dusty) 
plasma physics show that this system exhibits a com- 
plicated behavior which depends on the behavior of its 
"subsystems" which are represented by electrons, ions, 
neutral gas, and charged dust particles. All of them have 
their own relaxation time to local equilibrium; therefore, 
interaction among them may lead to a delay of relaxation 
to local equilibrium in plasma. Moreover, in the elec- 
tronic subsystem of plasma, local equilibrium does not 
exist, that stimulates development of theories beyond lo- 
cal equilibrium [o4|. Interaction of different subsystems 
in complex (dusty) plasmas with missing local thermody- 
namic equilibrium in the electronic subsystem, makes the 
description of observed experimental data of this object 
rather complicated. 

Experimental results of Morfill et al. |65| from plasma 
observations exhibit unusual behavior from weak colli- 
sionless interaction of gases to fluid flow with further 
possible crystallization of plasma. These results are de- 
scribed by means of molecular dynamic simulations |65| . 
The field approach, also, seems to be also applicable due 
to the fact that during transitions in plasma, the char- 
acteristic size of patterns is on the mesoscopic or even 
macroscopic scale. The field approach to a heat- and 
electronically-conducting fluid has been demonstrated in 
ionized gases |5l| using equations of generalized type of 
Eqs. O-Jl)^ 

VII. CONCLUSIONS 

The diffuse-interface model for rapid phase transfor- 
mation in metastable binary systems has been presented. 
To describe the steep but smooth change of phases within 
the width of diffuse interface, we use the formalism of the 
phase-field model. 

Rapid phase transformations may proceed under lo- 
cal nonequilibrium conditions. In our phenomenological 
macroscopic description, we extend the classic set of in- 
dependent thermodynamic variables by inclusion of dis- 



sipative fluxes as additional basic variables. The evolu- 
tion of the fluxes is characterized by their own dynamics 
with relaxation times r summarized in Table [i] Thus, 
the extended set Q of variables allows one to describe 
phase transformations with finite interface velocity that 
is comparable or even higher than 1/t, where I is the 
mean- free-path of particles (atoms). 

The evolution equations for the hyperbolic phase-field 
model with dissipation are derived from an entropy func- 
tional (|14f> based on the extended set Q of independent 
thermodynamic variables. This model yields a definite 
positive entropy production ((27(1 in consistency with the 
second law of thermodynamics. 

Generalization of the model has been done by introduc- 
ing memory functions and using a variational principle. 
As a result, the consistency of the macroscopic approach 
with the microscopic fluctuation-dissipation theorem has 
been found for the phase-field with memory [Eqs. 134(1 - 
1(37(1 ] : and nonlinear evolution equations [Eqs. 1(44(1 - 1(46(1 ] 
are derived from the variational principle 1(41(1 . 

The derived equations for the evolution of diffuse inter- 
face were correlated with existing models of nonequilib- 
rium transport processes and for systems under phase 
transformation. Particularly, we compare our deriva- 
tion to models of superconductivity, phase separation, 
viscoelastic or electronically-conducting fluids, interface 
motion by mean curvature, rapidly solidifying systems, 
and reaction-diffusion systems. 
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